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Abstract 

We develop a multidimensional coupled channel method suitable for studying the interplay of 
bound state resonance and phonon assisted scattering of inert gas atoms from solid surfaces in 
one, two and three dimensions. This enables us to get insight into the features that depend on 
the dimensionality of inelastic resonant processes typically encountered in low energy He atom 
scattering from surfaces, in general, and to elaborate on the observability of recently conjectured 
near threshold resonances in scattering from Einstein phonons, in particular. 
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A unified treatment of inelastic, resonant and diffractive scattering of inert gas atoms 
from solid surfaces has remained a challenging task since the very beginning of development 
of quantum atom-surface scattering theory in the 1930si. Already at that stage it was 
noted that a class of processes involving emission and absorption of phonons combined with 
the intermediate projectile propagation through bound states of the static atom-surface 
potential may take place and thereby modify the scattering amplitudes. However, in the 
works quoted in Ref. Q the descriptions of inelastic scattering and accommodation were 
limited to one-phonon processes that were treated essentially within the one-dimensional 
(ID) framework, i.e. only the component of the projectile momentum normal to the surface 
was assumed to change upon phonon emission or absorption. The three dimensional (3D) 
treatment was implemented only in the description of diffraction and selective adsorption 
processes. 

Further developments of the early quantum atom-surface scattering theories proceeded 
mainly in three directions: (i) full 3D unitary treatment of diffraction of atoms from 
surfaces^, (ii) full 3D unitary treatment of one-phonon inelastic scattering 3 and its extension 
by inclusion of diffraction^, and (iii) quantum description of sticking and accommodation 
processes^, dominantly within the one-phonon theories^, with extension to the selective 
adsorption processes^. A more complex problem of multiphonon atom-surface scattering, 
important for instance for the determination of ubiquitous Debye- Waller factor (DWF) and 
pronounced satellite structures in the scattering spectra, was tackled a bit later and mainly 
within the semiclassical theories 9,10 . Only recently a tractable fully quantum scattering spec- 
trum formalism for treating multiphonon atom-surface scattering has been developed^ 1 and 
applied to the interpretation of inelastic He atom scattering (HAS) spectra for the various 
systems in which the diffractive and inelastic spectral features could be easily disentan- 
gled. Although the scattering spectrum formalism offers exact solutions also in the case of 
a strong interplay between multiphonon inelastic and elastic resonant processes (i.e. diffrac- 
tion and selective adsorption) its applicability in this situation is hindered by computational 
difficulties. 

Diffractive scattering of atoms from ordered solid surfaces can be efficiently treated in 
the absence of inelastic processes by the method of coupled channel (CC) equation a 12 i 13 i 14 i 15 . 
Application of the CC method to the calculations of diffraction spectra for He atoms scat- 
tered from Xe monolayer adsorbed on graphit o 1 ^ 1 ^ demonstrated the potentiality of the 
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method in the treatment of surface scattering. Partial inclusion of the effect of substrate 
phonons through the Debye- Waller attenuations of diffraction peak intensities 1 ^ then posed 
a question of the extension of the CC approach so as to enable the assessment of resonant 
and inelastic effects and their interplay on an equal footing. This is possible in principle 
but difficulties arise in the implementation of such a demanding numerical algorithm to 3D 
scattering in which the phonon configuration space is multi-dimensional 13 . On the other 
hand, the treatment of ID inelastic scattering by Einstein phonons (colinear motion of the 
projectile and surface oscillator) by the CC method poses no difficulties 14 and it was first 
applied in atom-surface scattering to investigate the Debye- Waller factors^ and the long 
standing problem of the importance of linear and nonlinear projectile-phonon coupling in 
inelastic scattering in the absence of bound states 2 ^. Extension of this analysis by the in- 
clusion of bound state resonance s 21 > 22 showed that the latter can indeed significantly affect 
the amplitudes of ID low energy scattering near the single-phonon excitation threshold, in 
a sequence of strongly correlated events of single- and multiphonon excitations and interme- 
diate projectile propagation through one or more bound states. These findings have posed 
a fundamental question as to whether such a practically exact unified ID treatment of res- 
onant and inelastic scattering can predict and interpret some new features in the 3D high 
resolution- low energy angular resolved HAS from surfaces 2 ^. 

The answer to the above posed question obviously requires extension of the CC method to 
description of inelastic scattering beyond ID. This is almost a formidable task if the phonon 
configuration space is multi-dimensional, i.e. if the phonon quantum numbers (wavevector 
and polarization) span a multi-dimensional Fock configuration space. To overcome this diffi- 
culty we introduce here a tractable restricted Fock space coupled channel (RFCC) algorithm 
for treating inelastic atom-surface scattering involving bound state resonances that applies 
equally well to 1-, 2- and 3D collisions. 

A prerequisite for accurate description of low-energy inelastic atom-surface scattering is 
a correct treatment of the momentum transfer that takes place between the projectile and 
substrate phonons. To implement this requirement in a numerical treatment we begin in a 
standard fashion and discretize the phonon wave vector space 1 ^-. We take N discrete values 
(Qi? Q2, Qn) of phonon wavevector Q that are distributed homogeneously over the first 
surface Brillouin zone. If one takes into account the states with j phonons excited, the 
dimension of the basis set describing the subset of these states grows as NK Owing to 
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this, a numerical treatment of higher excited phonon states j > 2 for D > 1, even with a 
sparsely discretized space of phonon wave vectors, is beyond the capabilities of the present 
day computers. Hence, in the present RFCC approach we consider only the phonon ground- 
state, the subset of N singly excited states, and the subset of N(N — l)/2 doubly excited 
states, i.e. our restricted Fock space of excited states contains (iV 2 + N + 2)/2 different 
phonon states that should suffice for describing typical regimes of HAS experiments carried 
out to reveal the characteristics of surface phonons in the single phonon excitation regime^ 
and beyond 24 ! 25 ' 26 . 

We start from the total wave function of the system 

(R, z, {u,}) = 4>n{z) exp(-zK n R) ® \n) (1) 

n 

where r = (R, z) is the radiusvector of the projectile (of mass m), with R lying in the surface 
plane, the coordinates {u^} describe the displacements of the target atoms of mass M, I is 
the lattice site index, K n is the projectile wavevector parallel to the surface, and |n) is a 
Fock state describing phonon excitations in the target. The functions <p n (z) satisfy the CC 
equations 

d 2 

— (f) n (z) = ^ W n , n >{z)(f) n >{z), (2) 

n' 

and obey the scattering boundary conditions 17 . Henceforth we shall restrict the coupling to 
z-polarized displacements only and assume normal projectile incidence to enable equivalent 
treatment of the scattering event in all dimensions. The CC matrix elements are given by 

W n ,Az) = ^V^(z)5 n>n , - ^ n ,n< + ^jrV*$(z). (3) 
The reduced channel energy, k 2 , is 

2 ™( Ei -E n )-Kl (4) 
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where E{ is the projectile incident energy, E n = J2 p =i n Q p ^(Qp) is the vibrational energy 
of the target in the state |n) with uq p denoting the occupation of Q p th phonon mode 
in the state |n), and u;(Q) is the phonon frequency. Following the common practice we 
take the projectile-target atom binary potentials in separable form^. v(r — r^) = Vd{z — 
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zi) exp [— (R — Hi) 2 /(2a 2 )}, which when pairwise summed over I yield for the average total 
static interaction potential: 



where a is the nearest neighbor distance in the simple Bravais lattice of the target surface. 
The matrix of the dynamic part of the interaction potential is given by 



For computational convenience we have assumed only linear coupling of the projectile to the 
target atom displacements as this gives the dominant contribution to inelastic scattering 2 ^ 
The above equations do not incorporate transitions involving the reciprocal lattice vectors 
G 7^ i.e. they are strictly applicable to the case of smooth surfaces with 2ttcx ^> a for 
which the diffraction probabilities are negligible. The resulting set of CC equations is then 
solved using standard technique o 1211 ^ 141 " 1 ^ 1 ^ 17 . 

The above outlined RFCC formalism is applied to the paradigmatic case of He atom 
scattering from a vibrating Xe monolayer adsorbed on a cold Cu(lll) surfac o 25 ! 26 ' 27 ' 28 . The 
manifestations of inelastic bound state resonance effects are expected to be most pronounced 
in this prototype system because (i) the corrugation of the static He-Xe/Cu(lll) potential 
gives rise to extremely weak diffraction peaks and hence the G-assisted transitions can be 
neglected, and (ii) among the three Xe- monolayer localized phonon modes the vertically 
polarized S-mode, which couples most strongly to the scattered He atom, exhibits no dis- 
persion practically over the entire surface Brillouin zone, i.e. 7tu;(Q) = hu = 2.7 meV. 
Hence, we can restrict our RFCC-analyses of the bound state resonance effects to the ex- 
perimentally relevant situation in which He atoms are scattered by the vertically polarized 
Einstein phonons that propagate along a smooth surface. 

We represent the binary potential Vr>(z — Zq) by a Morse form, and adjust its well depth 
to obtain the same static component V stat (z) = Vq (exp(—2z/b) — 2exp(—z/b)) of the inter- 
action potential for all D. This gives rise to the transitions between the same energy levels 







5 




20 40 60 80 100 

E. [meV] 



FIG. 1: Comparison of two-phonon RFCC and multiphonon EBA results for the Debye- Waller 
factor Pq in 2D He atom-surface scattering as a function of normal incident energy E{. Model 
parameters are: Huq = 2.7 meV, a = a = 2 A. 



so that the corresponding probabilities can be directly compared. The Morse parameters of 
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and 



22. i-e. V Q = 6.6 meV and b = 0.82 A, yielding 



V stat (z) are the same as in Refs. 
eo = —4.54, ei = —1.67, and €2 = —0.14 meV. This choice facilitates comparison of the 
earlier calculated ID results with the present ones. Note, however, that our ID results can 



22j because in the present study the 



not be identical to those obtained in Refs. |21| and 
phonon Fock space is restricted and the model does not incorporate nonlinear coupling of 
the projectile to surface vibrations. As these restrictions are not easily removed in higher 
dimensions, we retain them here also in the ID calculation to enable consistent comparison 
of the results obtained for different dimensions. 

The reliability of the present two phonon RFCC approach is tested by comparing the 
elastic scattering probabilities Pq (i-e- DWF's) calculated by this method with the ones 
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FIG. 2: Comparison of the RFCC Debye- Waller factors Po arising from interplay of Einstein 
phonon exchange and bound state resonance effects in ID, 2D and 3D, and the Po obtained in 
3D EBA (dashed line) for multiphonon He atom scattering. The resonances appear around E^ = 
e± + 2fiujo = 3.83 meV and Ei = e<i + 2T%ujq = 5.3 meV. 

obtained from the multiphonon exponentiated Born approximation (EBA) 11 ' 19 that is exact 
in the low and high incident energy limits for all D and produces accurate off-resonance 
results also at intermediate energie a 19 i 20 i 21 i 22 . Figure ^ compares Po obtained from the RFCC 
and EBA calculations for the intermediate case D = 2. It should be noted that numerical 
requirements on these two types of calculations are quite different. The EBA formalism 
requires calculation of Q-resolved single-phonon excitation probabilities and their integration 
over the first Brillouin zone^. On the other hand, the RFCC calculations involve numerical 
propagation of large matrices [whose number of rows and columns is (A^ 2 + N + 2)/2] along 
the z-coordinate, and are thus much slower from the computational aspect (for the 3D case 
illustrated in Fig. El they are slower by a factor ~ 10 5 ). The displayed EBA and RFCC 
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results significantly differ only in the low energy region where Pq is affected by the resonant 
projectile propagation through the bound states by emission and reabsorption of phonons. 
Inset shows the region of high incident energies in which the validity of the 2D RFCC 
calculation breaks down due to the neglect of Fock states with more than two excitation 
quanta. In the EBA, the average number of phonons n excited in a scattering event is related 
to Pq through P = exp(— n) 11 . Thus, the displayed two phonon RFCC results are expected 
to be reliable in the scattering regime in which P$ > exp(— 2) = 0.135, i.e. above the dotted 
line in the inset. The EBA predicts that this regime is reached for E{ ~150 meV, which 
is precisely where the RFCC results begin to exhibit pronounced oscillatory behavior and 
deviation from the corresponding asymptotically exact multiphonon EBA values. Hence, an 
excellent agreement between the results for off-resonant scattering in the range < E± < 150 
meV obtained from two completely different algorithms confirms the validity of the present 
two-phonon RFCC approach to study the interplay of multiphonon and resonant scattering 
at low Ei. This is crucial because the resonant scattering effects predicted in the RFCC 
approach cannot be encompassed by the simplest version of the EBA even in the ID cas o 21122 . 

The thus established validity of the two phonon RFCC for treating the interplay of mul- 
tiphonon and resonant effects in the experimentally relevant case of He atom scattering 
from Einstein phonon s 25 ! 26 enables the assessment of the effect of dimensionality on these 
processes. Figure El displays a comparison of the elastic scattering probabilities Pq for He 
atoms incident normal to the "surface" of a ID, 2D and 3D system. These results demon- 
strate that the inelasticity of scattering and the bound state resonance effects are reduced 
as D increases. In the ID case of colinear scattering the on-the-energy and momentum- shell 
requirements allow only the inelastic transitions for which Q = 0, and hence they are all 
governed by the maximum values of V dyn (cf. Eq. EJ). On the other hand, with the increase 
of D the on-shell scattering intensities are redistributed over the parts of momentum space 
with nonzero Q in which the magnitude of the matrix elements (jBJ) is reduced. As a result, 
as D increases the total contribution to inelastic and resonant scattering intensities dimin- 
ishes. Hence, although the signature of inelastic bound-state resonances persists in higher 
dimensions, it is already too weak in 3D to be experimentally observed for smooth surfaces 
on which the G-assisted transitions are negligible. 

In summary, we have developed a restricted Fock space coupled channels method suitable 
for a unified treatment of the interplay between multiphonon excitation and bound state 
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resonance effects in atom-surface scattering. Application of the method to a prototype 
system of He atom scattering from Einstein phonons on smooth surfaces of adsorbed Xe 
monolayers shows that, in contrast to the earlier conjectures based on ID models, such 
effects are only weakly discernible even under the most favorable experimental conditions 
for their observation^. 
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